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\^ ■ ABSTRACT 
ON ■ 

GN ' The charm physics program accessible to a Tau-Charm Factory is sum- 

marized. Semileptonic, leptonic, hadronic, rare and forbidden charm de- 
cay modes are discussed, as well as D° — D° mixing and CP violation 

/-»s ! in the charm sector. The theoretical expectations in the Standard Model 

and beyond as well as the experimental capabilities of a Tau-Charm Fac- 

ON ! tory are examined for each of the above processes. 

1 . Overview 

> . 

One of the outstanding problems in particle physics is the origin of the fermion 
£T) ■ mass and mixing spectrum. At present the best approach in addressing this question 

is to study the properties of all heavy fermions in detail. While investigations of the K 
and B systems have and will continue to play a central role in the quest to understand 
flavor physics, in-depth examinations of the charm-quark sector have yet to be per- 
& ' formed, leaving a gap in our knowledge. Since charm is the only heavy charged +2/3 

quark presently accessible to experiment, it provides the sole window of opportunity to 
examine flavor physics in this sector. In addition, charm allows a complimentary probe 
of Standard Model (SM) physics (and beyond) to that attainable from the down-quark 
j> | sector. 

| Detailed measurements of heavy quark systems are best realized at high precision, 

high luminosity machines. Here we discuss the charm physics potential at a high lumi- 
nosity (C = 10 33 cm -2 sec -1 ) e + e~ collider operating at the cc and rf threshold region, 
i.e., a Tau-Charm Factory (rcF). Such a machine would facilitate an in-depth analysis 
of the charm-quark sector (as well as r-leptoni and charmoniumi studies) without the 
contamination of 6-quark production. The high luminosity is absolutely essential in 
order to accomplish the physics goals of the machine. In some cases it is necessary in 
order to measure a reaction at the high level of precision that is desired, while in other 
cases, high luminosity is required to reach the very small transition rates. 
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y/s GeV 


Particle 


Cross Section (nb) 


Pairs Produced (xlO 8 ) 


3.77 


D°D° 


5.8 ±0.8 


1.0 


3.77 


D+D- 


4.2 ±0.7 


0.8 


4.03 


D S D S 


0.7 ±0.2 


0.24 


4.14 




0.9 ±0.2 


0.32 



Table 1. Charm production cross sections and rates. 



The exclusive nature of charm production at a rcF automatically supplies the 
production kinematics and low combinatoric backgrounds that are essential to back- 
ground rejection. The charm production cross section in this energy region is large and 
well measured, as can be seen in Table 1 from Ref. [3]. Although the largest production 
rates occur at threshold, the proposed 5-Factories and a high luminosity Z Factory 
could produce a comparable (but slightly smaller) number of charmed particles. The 
advantage of the threshold region at a rcF is that it provides the capability to control 
backgrounds and systematic errors. Since the beam energy can be tuned to lie just 
below or above production threshold, physics backgrounds can be directly measured, 
instead of estimated via Monte Carlo simulations. Because heavy flavors are pair pro- 
duced at threshold, the observation of the decay of one particle, cleanly tags its partner. 
This yields a sample! of 10 6 — 10 7 singly tagged D-mesons per year, a size which is 
crucial in order to obtain highly precise measurements of branching fractions and to 
search for rare processes. 

The charm physics program at a rcF is strong and diverse. Charm hadrons pos- 
sess a rich variety of weak decays: Cabibbo allowed, single Cabibbo suppressed, double 
Cabibbo suppressed, leptonic, semileptonic, and rare higher-order decays. The scale of 
the interactions is in the middle of the regime where perturbative and non-perturbative 
QCD effects are cleanly separable. Hence the charm system provides an excellent labo- 
ratory to test our understanding of QCD, the dynamics of heavy quark decay, and the 
structure of the hadronic weak current. The knowledge gained from a comprehensive 
study would supply the needed benchmarks for lattice QCD and result in a more reli- 
able extrapolation to the B system. Fundamental parameters in the SM, such as the 
values of the Cabbibo-Kobayashi-Maskawa (CKM) mixing matrix elements involving 
charm-quark decay, would also be more precisely determined. The charm system also 
offers the opportunity to study CP violation in the up-quark sector and to search for 
new sources of CP violation. Loop induced processes are very sensitive to new particles 
which may participate in the loop and hence may provide a signature for physics be- 
yond the SM. Rare charm interactions compliment searches in the down-quark sector 
because the couplings to new particles may be flavor dependent, either through mass 
dependent couplings or through mixing angles. 

We now turn to a discussion of semileptonic, leptonic, hadronic, rare and forbid- 
den decays of the charm quark, as well as D° — D° mixing and CP violation in the 
charm sector. In each case we (i) review the motivation for studying these processes, 
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(ii) summarize the present status of theoretical calculations in the SM and beyond, and 

(iii) examine the experimental sensitivities available at a tau-charm factory. Further 
details can be found in Refs. [4,5]. 

2. Semileptonic Decays of Charmed Hadrons 

Semileptonic decays occupy a fundamental place in the experimental study of 
charm decays and in attempts to understand the underlying dynamics associated with 
charmed hadrons. They also play a critical role in determining the elements of the CKM 
mixing matrix. However, the extraction of the values for the CKM elements requires 
knowledge of the form factors which parameterize the matrix elements of the hadronic 
weak current. Theoretical calculation of these form factors is a thorny problem, particu- 
larly in the case of charm decays as they occur at an energy scale which lies between the 
regimes where chiral perturbation theory and heavy quark mass expansions (which are 
used to describe K and B decays, respectively) are absolutely valid. With the excep- 
tion of lattice calculations, the form factors cannot be calculated from first principles, 
and thus model dependent estimates, usually combined with approximate symmetry 
properties, must be employed. It is therefore crucial to measure the shape of the form 
factors in order to test the theoretical models. This would enhance our understanding 
of QCD, and enable scaling of the predictions to the B sector with better accuracy, as 
well as provide an improved determination of the CKM matrix elements in the charm 
sector. 

Before turning to a discussion of the form factors, we first note that the CKM 
elements governing semileptonic charm decay, V c d and V cs , are surprisingly poorly deter- 
mined at present. The best determination of V c d is from charm production in neutrino 
scattering off valence (i-quarks, yielding the value! \ V c d\ = 0.204 ±0.017. Extracting V cs 
from this method is plagued by uncertainties from the estimates of the strange quark 
parton density and hence the semileptonic decay D + — > Ke + v is used. Comparison 
of the data with theory and using conservative assumptions in calculating the form 
factors yields! \V CS \ = 1.01 ± 0.18. The imposition of CKM unitarity for a three gen- 
eration SM constrains! the values of these elements further to \V c d\ = 0.221 ± 0.003 
and | | = 0.9745 ± 0.0007. We see that these CKM elements are very well known 
if one assumes unitarity, but that the direct experimental measurements are at the 
10 — 20% level. There are many theories beyond the SM which preclude a three gen- 
eration unitary CKM matrix and they should be tested with a precision in the charm 
sector that is equivalent to that of planned experiments in the B system. For example, 
the comparison of the exclusive semileptonic decays D — > i\{iv and D — > Kfiu, where 
theoretical uncertainties mostly cancel in taking the ratio, would allow a determina- 
tion of | | / 1 | to ~ 1%, which is comparable to the present level of precision on the 
Cabibbo angle 9 C . 

Semileptonic rates are dominated by modes containing only one hadron in the 
final state. These are K or K* for Cabibbo allowed decay, 7r or p for Cabibbo suppressed 
D decay, and 77, 77', for Cabibbo suppressed D s decay. For the decay into a pseudoscalar 
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meson, D — > P£v } the hadronic matrix element can be written as 

(P(k)\q lfl c\D(p)} = f + (q 2 )(p + k), + f4q 2 )(p - , (1) 

where q 2 = (p — k) 2 is the momentum transfer and q = s,d for P = K,n, respectively. 
If the lepton mass is neglected, only / + contributes to the q 2 distribution with 



dT _ G 2 F \V cs(d) \ 2 k 3 
dq 2 24vr 3 



cs(d) | K r , / 2\ 1 2 



\U(q 2 W- (2) 



We note that an accurate comparison of the rate for D — > Kev to that for D — > K\iv 
would determine f-{q 2 ) and provide the best direct measurement of V cs . 

Theoretical uncertainties are introduced by assuming a particular q 2 dependence 
of the form factor. It is usually parameterized using vector meson dominance with 

Utf) = : /+ 2 ( ° } . 2 , (3) 
1 — q z /m* z 

where m* is the pole mass of the appropriate channel. This method allows one to 
fit and determine the pole mass. In heavy quark effective theory (HQET), all decay 
modes are described in terms of a single form factor assuming an exponential q 2 depen- 
dence, f+(q 2 ) = /+(0) exp(ag 2 ). This yields values of /+(0) which are approximately 
5% smaller than those obtained in the single pole ansatz. However, HQET is less ad- 
vantageous in D decays as the symmetry breaking corrections as expected to be large. 
Various theoretical predictions^ (from quark model calculations, QCD sum rules, and 
lattice QCD) for the form factor in D — > Klv, as well as the average experimental 
result ,i deluding a recent determination by CLEoi from a measurement of the f dis- 
tribution, is summarized in Table 2. In the case of the monopole distribution, the pole 
mass is then measured to be m* = (2.00 ± 0.12 ± 0.18) GeV, which is consistent with 
the D* mass. For HQET, the fit value of a is a = (0.29 ± 0.04 ± 0.06) GeV" 2 , which 
should be compared with the predicted value! of a = 0.21 GeV -2 . 

The various quark model predictions for the q 2 dependence are in rough agreement 
with each other as they all incorporate a recoil energy dependence that reflects a fall-off 
with a typical hadronic scale of ~ 1 GeV. Large disagreements in the predictions will 
only occur either at low energies, if there are poles close to the physical region, or at 
recoil energies greater than 1 GeV. Therefore charm meson semileptonic decays do not 
fully probe the different aspects of the models given the limited kinematic range. Very 
precise measurements will then be needed in order to discriminate amongst the models. 

For semileptonic decay to a vector meson V, the hadronic matrix element can be 
parameterized in terms of four form factors, 

(V{k)\q llx c\D(p)) = 2V{ / ] e^e'Vfc* , (4) 

m D + m v 
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Reference 


f fr\\ 
/+l u ) 


1/(0) 






Exp. Average 


n _i_ n no _l n no 
[).(& ± U.Uz ± U.Uz 


ii i n o 
1.1 rh U.z 


U.00 ± U.U4 


U.4U ± U.U8 


1 heory 










TO r^WT 

1SGW 


m on 

(J. 82 


l.i 


0.8 


0.8 


BSW 


0.76 


1.3 


0.88 


1.2 


AW 


0.7 


1.5 


0.8 


0.6 


BKS 


0.90 ±0.08 ±0.21 


1.4 ±0.5 ±0.5 


0.8 ±0.1 ±0.3 


0.6 ±0.1 ±0.2 


LMMS 


0.63 ±0.08 


0.9 ±0.1 


0.53 ±0.03 


0.2 ±0.2 


BBD 


0.60 


1.1 


0.5 


0.6 



Table 2. Measured and predicted form factors in D — > and -D — > K*lv. 



Reference 


T(D - 


-> K*iu)/T(D -> 


Experiment 


0.58 ±0.06 


ISGW 




1.14 


BSW 




0.87 


AW 




1.34 


BBD 




0.6 ±0.4 



Table 3. Theoretical predictions and experimental average for the ratio T(D — > K*£v)/T(D — > 



A (a 2 

[m D + m v )A l {q 2 )e* 2 y e* ■ p(p ± fc)„ 

^ m c + my 

±^A 3 (g 2 )6*.p(p-^ 



where e represents the polarization vector of V. The contribution from A^q 2 ) van- 
ishes in the limit of massless leptons. Here, branching fraction measurements alone do 
not allow the determination of the form factors, and angular distributions must also 
be measured in order to extract all three. A summary of existing form factor mea- 
surements! and various theoretical predictions! is displayed in Table 2 for comparison 
purposes. We see that most calculations agree, within errors, with the experimental 
determination of V^(0), but fail to predict Ai(0) and A 2 (0). 

Theoretical predictions! for the ratio of vector to pseudoscalar decay modes are 
compared to the experimental averaged result! in Table 3. We see that the QCD sum 
rule approach (BBD) obtains the best agreement with experiment, but also has large 
errors. The agreement with the /+(0) measurements seen in Table 2 points to the K* 
mode as the problem, which is also in accord with our conclusions on the axial vector 
form factors in Table 2. Since we seem to lack a clear theoretical understanding of these 
decays, accurate experimental information in several modes is vital. 

There are two recent experimental measurements of the vector and axial-vector 
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Reference 


wo)Mi(o) 


A 2 (0)/A 1 (0) 


Exp. Average 


2.0 ±0.7 


1.8 ±0.5 


SI 
BKS 
LMMS 


1.85 

2.00 ±0.19 ±0.23 
1.65 ±0.2 


1.21 

0.78 ±0.08 ±0.15 
0.33 ±0.36 



Table 4. Measured and calculated form factors in D s — ► </>£z/. 



form factors in D s decays. They constitute a further opportunity to examine the dis- 
crepancies between theory and experiment which appear to exist in D decay. Table 4 
summarizes the measured! and predicted! form factors in the decay D s — > <f>i,v. We see 
that theory fails to provide the correct value for the ratio A 2 (0)/Al(0), just as in the 
K* mode. It would be of interest to investigate the process D s — > (r/ + r/)£z/ to see if 
the disagreement between theory and experiment persists only in the vector mode as 
in the case of D decays. 

Charmed baryon semileptonic decays provide yet another opportunity to test the- 
oretical calculations. For example, a measurement of A+ — > Aiu, which is expected to 
dominate the A c semileptonic decays, would yield information on heavy quark predic- 
tions that relate these decays to D — > X s iv. Measurement of the A c and S c branching 
factions would supply essential input to the analysis of the charm content in B decays. 

A study of inclusive semileptonic decays would provide an additional test of 
HQET, as well as allowing for the extraction of a value for the charm quark mass. The 
results of a recent HQET analysis^ from present data on charm semileptonic decays 
yield m c = 1.57 ± 0.03 ± 0.05 GeV, where the errors represent the experimental and 
theoretical uncertainties, respectively. Another important question to be addressed is 
whether the sum of the exclusive rates saturate the inclusive rate. 

With adequate statistics, the exclusive semileptonic rates as well as the q 2 de- 
pendences and relative normalization of the form factors can be measured at a rcF. 
The ability to make efficient measurements over the full kinematic range is essential. 
The expected number! of detected exclusive modes per year at a rcF lies in the range 
10 4 — 10 5 for the various channels in D°, D + decay and 10 3 — 10 4 for D s . Figure 1 from 
Ref. [10] displays the use of tagging, kinematics, and calorimetry to isolate Cabbibo 
allowed and suppressed decays by the method of reconstructing the missing mass. This 
analysis shows that the events are expected to be essentially background free. A rcF 
could provide an exhaustive and accurate data base, with which the various theoreti- 
cal models could be compared, and allows for the reliable determination of the CKM 
elements governing charm-quark decay. 
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Meson 




T + U T 


D + 


3.52 x 10" 4 
4.21 x 10~ 3 


9.34 x 10" 4 
4.11 x 10~ 2 



Table 5. SM branching fractions for the leptonic decay modes, assuming fo 
f D „ = 230 MeV. 



200 MeV and 



3. Leptonic Decays of Charmed Mesons 

The SM transition rate for the purely leptonic decay of a pseudoscalar charm 

riD+^e+^^^fnjV^mn^mjll-^) , (5) 



meson is 



m 



d 



with q = d,s and fo ( q) is the weak decay constant defined as usual by 



(0| W5e|£>( 9 )(p)) = ifDr q) P» 



(6) 



where f„ = 131 MeV in this normalization. The resulting branching fractions are small 
due to the helicity suppression and are listed in Table 5 using the central values of the 
CKM parameters given in Ref. [6] and assuming fr> = 200 MeV and fo B = 230 MeV. 

Assuming that the CKM matrix elements are well-known, the leptonic decays can 
provide important information on the value of the pseudoscalar decay constants. Pre- 
cise measurements of these constants are essential for the study of D° — D° mixing, CP 
violation in the charm sector, and non-leptonic decays. They would also test the accu- 
racy of QCD calculational techniques and provide a means for more precise predictions 
when extrapolating to the B sector. The existing upper limit for is fo < 290 MeV, 
and is derived from the 90% C.L. bound§ B(D + pL+v^) < 7.2 x 10" 4 . CLEO has 
recently observedS the process D* + — > Df'j — > fiwy by examining the mass difference 
5M = M^ U1 — and have obtained 



T(D+ -> 07T+) 



0.235 ±0.045 ±0.063 



(7) 



Using F(D+ -> 07T+) = 3.7 ± 1.2% they find f Da = 337 ± 34 ± 45 ± 54 MeV where 
the last error reflects theuncertainty in the 07r + branching fraction. An emulsion 
experiment has measuredli2l f Ds = 232 ± 45 ± 52 MeV. The BES Collaboration has 
also recently reporte d0 the observation of three candidate events in e + e — > DfD s 
with the subsequent decay D s — > fiu yielding fjj a = 434^^3 t| 3 MeV. The errors are 
expected to improve once more statistics are obtained. 

A variety of theoretical techniques have been employed to estimate the value of 
fr> and fo s . Lattice QCD studiesS calculate these quantities in the quenched approxi- 
mation through a procedure that interpolates between the Wilson fermion scheme and 
the static approximation. The non-relativistic quark model is used to relate the decay 
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Decay Constant 


Lattice 


Quark Model 


Sum Rule 


fn 


208(9) ± 35 ± 12 


207 ± 60 


170 - 235 


/d s 


230(7) ± 30 ± 18 


259 ± 74 


204 - 270 




1.10(2) ±0.02 ±0.02 ±0.03 


1.25 


1.21 ±0.06 



Table 6. Theoretical estimates of the weak decay constants in units of MeV, (taking m c 
1.3 GeV in the sum rule approach). 



constant to the meson wave function at the origin, /m = x/^/Mm ^(O)!, which is 
then inferred from isospin mass splitting of heavy mesons.Ej Another approach uses 
QCD sum rules.0 For each of these calculational methods, the resulting ranges for 
the values of the pseudoscalar decay constants are presented in Table 6. Given the 
large errors, the results of these three approaches are roughly consistent. The theo- 
retical uncertainties associated with the ratio /d«//d are much smaller, as this ratio 
should deviate from unity only in the presence of broken SU(3) flavor symmetry. The 
magnitude of such flavor violating effects can be determined by the measurement of 
B(Df — > /i + ^ M ) / B(D + — > /U + z/ M ). Obtaining precise values for these decay constants, 
both theoretically and experimentally, would provide valuable input for a wide range 
of phenomenological applications. 

Non-SM contributions may affect the purely leptonic decays. Signatures for new 
physics include the measurement of non-SM values for the absolute branching ratios, 
or the observation of a deviation from the SM prediction 



B(D+ S) - p+vj mj (l - ml/ml J 



B(D 



mi 



m^/m 2 D 



(8) 



This ratio is sensitive to violations of \i — r universality. 

As an example, we consider the case where the SM Higgs sector is enlarged by 
an additional Higgs doublet. These models generate important contributions!!! to the 
decay B — > tv t and it is instructive to examine their effects in the charm sector. 
Two such models, which naturally avoid tree-level flavor changing neutral currents, are 
Model I, where one doublet (02) generates masses for all fermions and the second (0i) 
decouples from the fermion sector, and Model II, where 02 gives mass to the up-type 
quarks, while the down-type quarks and charged leptons receive their mass from 0i. 
Each doublet receives a vacuum expectation value v iy subject to the constraint that 
v i + v 2 = v sm- The charged Higgs boson present in these models will mediate the 
leptonic decay through an effective four-Fermi interaction, similar to that of the SM 
H^-boson. The interactions with the fermion sector are governed by the Lagrangian 



C 
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2V2M W 



-H ± [V ij Tn u .A u Ui(l - j 5 )dj + Vijm d] A d Ui{\ + j 5 )dj 
miAtVi(l + 7 5 )£] ± h.c. , 



(9) 
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with A u = cot (3 in both models and A d = Ag = — cot /3(tan/3) in Model I(II) , where 
tan/? = v<ijv\. In Models I and II, we obtain the result 

B ( D +^l+ Ue ) = B SM (l + ^f-) , (10) 

V m H±J 

where in Model II the decay receives an additional modification 



B(D+ - t+v t ) = B 



SM 



2 / m ^2 



1 + —p- 1 - tan' (3— 



(11) 



In this case, we see that the effect of the H ± exchange is independent of the leptonic final 
state and the above prediction for the ratio in Eq. (8) is unchanged. This is because the 
H contribution is proportional to the charged lepton mass, which is then a common 
factor with the SM helicity suppressed term. However, the absolute branching fractions 
can be modified; this effect is negligible in the decay D + — ► £ + U£, but could be of order 
a few percent in Df decay if tan /3 is very large. 

The detection of D + ,Df — ► /i + z/ decays is straightforward with single tagged 
event samples at a rcF allowing the measurement of the weak decay constants at the 
1% levelHl The t + u modes can also be observed via the subsequent decays r — ► zu, Ivv 
by relying on tagging, hermiticity, and background measurements. The tagged event 
samples are used to provide a constrained fit for the neutrino mass as well as to suppress 
backgrounds. Figure 2 from Ref. [19] shows that the missing mass spectrum cleanly 
separates the signal from background in all cases. 

4. Nonleptonic Decays of Charm Hadrons 

Nonleptonic charm decays provide another opportunity to test our theoretical 
understanding of heavy flavor decays. They occur in a regime where non-perturbative 
effects are more (less) important than in B (K) decays. The thorough investigation of 
the rich variety of decay processes available would unravel these QCD effects and would 
build a foundation for a more reliable extrapolation to the underlying dynamics in the 
B system. Several types of processes mediate nonleptonic charm decays, e.g., external 
spectator, color suppressed or internal spectator, weak annihilation, exchange, pen- 
guin, and mixing. Disentangling the underlying quark processes and determining their 
relative strengths is difficult and will require an increase of 2 — 3 orders of magnitude 
over current statistics, but will generate a detailed understanding of strong interaction 
effects. 

The absolute branching fractions of hadronic charm decays need to be measured 
with better precision. The present sensitivity of measurements of D° and D + decay is 
at the 10 — 20% level, while the situation for D s is even more uncertain as the absolute 
scale of its branching fractions is unknown. Most D s decays are measured relative to 
the 07r + mode, which is determined from the ratio T(D^ — > 0e + z/)/r(Z)+ — > 07r + ) 
combined with estimates of the cross section, theoretical predictions relating D + 
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and Df decays, and fragmentation assumptions. A direct determination of the <f)ir + 
branching fraction could thus shift the values of the remaining branching fractions. 
Very few of the rates for charm baryon decay (particularly for H c and Q c ) have been 
measured. Accurate data on these modes would provide invaluable assistance in the 
determination of the charm content in B decays, CP violating asymmetries, and in 
fragmentation studies. 

The rcF is well suited! for these measurements. The increase in statistics and the 
efficient single and double tagging techniques allows the determination of the absolute 
branching fractions for D° and D + at the 1% level, limited only by systematics. Ap- 
proximately 5 x 10 5 single D s tags can be reconstructed in one year at 10 33 cm -2 sec -1 
using the <f)ir + , S*tt + , K°K + and K*°K + decay modes, while the number of double 
tags from pairing these channels is 4950 x i?(</>7r + ). This yields a ~ 3% sensitivity to 
the absolute scale of D s decays. Similar absolute measurements of charm baryon decays 
could also be performed. 

D mesons are sufficiently heavy for many body final states to play a promi- 
nent role in their nonleptonic modes, however, D° and D + decays presently appear 
to be dominated by two body modes. In this case, most of the theoretical work is 
based on factorization^ or QCD sum ruleS approaches. The latter technique yields 
approximately model independent predictions for the amplitudes of four decay types 
D ,D S — > PP, PV, where P(V) is a pseudoscalar (vector) meson built of light quarks, 
but the numerical results depend on the assumption of SU(3) flavor symmetry. The 
former method is based on PP, VP, VV, and AP final states (where A represents an 
axial-vector meson) and assumes, (i) an amplitude which is expressed as the product 
of two current matrix elements, i.e., 

(M 1 M 2 \J fM J fl \D) ~ (Mil J"\0) ■ (M 2 \J^\D) , (12) 

(ii) two fit parameters, a\ and a 2 , which occur in the nonleptonic weak Hamiltonian 

H eff = JfF . [ ai (ud')(s'c) + a 2 {sd')(uc)\ : , (13) 

where the colons denote Wick ordering and d' , s' are Cabbibo rotated, (iii) model 
dependencies for the hadronic wavef unctions, and (iv) that the contributions from 
penguin amplitudes and final state interactions are negligible. This model currently 
provides a reasonable description for D° and D + decay modes, but may not hold up 
once more accurate data is obtained. The systematic study and comparison of D and B 
decays will help clarify some of the unresolved issues and inherent model dependencies 
associated with factorization. 

The situation with D s decays is not so well described. The anticipated strength of 
the weak annihilation process has yet to be confirmed and it is possible that the pattern 
of D s decays may not mimic those of D° and D + . It has been suggestedlH that the 
D s may have an enhanced non-resonant width into multi-particle final states. In this 



10 



case, the constraints of production at threshold help to reconstruct these complicated 
many body modes. An reconstruction example from Ref. [3] at y/s = 4.028 GeV for 
D s — ► ?77r + 7r + 7r~ with the subsequent decays r\ — > 7r + 7r _ 7r° and 7r° — > 77 is depicted in 
Fig. 3. 

Doubly suppressed Cabibbo decays, having branching fractions of order tan 4 9 C 
relative to Cabbibo allowed decays, also test our understanding of the hadronic weak 
current. These modes are best measured by tagging in D + decays, since there is no 
confusion with the possible mixing component which may be present in D° modes. A 
handful of events in D + — > K + K + K~ , <pK + and D° — > K + n~ have been observed.!^ It 
is expected! that a few hundred such events can be reconstructed at the rcF and that 
a signal sensitivity at the ~ 10% level can be achieved. 

5. Rare and Forbidden Decays of Charm Mesons 

Flavor changing neutral current (FCNC) decays only occur at the loop level in the 
SM. Due to the effectiveness of the GIM mechanism and the small masses of the quarks 
which participate inside the loops, short distance SM contributions to rare charm decays 
are very small. Most reactions are thus dominated by long distance effects which are 
difficult to reliably calculate. However, a recent investigation! of such effects indicates 
that there is a window for the clean observation of new physics in some interactions. 
In fact, it is precisely because the SM FCNC rates are so small, that charm provides 
an important opportunity to discover new effects, and offers a detailed test of the SM 
in the up-quark sector. 

FCNC decays of the D meson include the processes D° — > £ + £~,77, and D — > 
X + £ + £~,X + uu,X + 7, with £ — e,fi. They proceed via electromagnetic or weak 
penguin diagrams as well as receiving contributions from box diagrams in some cases. 
The calculation of the SM short distance rates for these processes is straightforward 
and the transition amplitudes and standard loop integrals, which are categorized in 
Ref. [22] for rare K decays, are easily converted to the D system. The loop integrals 
relevant for D° — > 77 may be found in Ref. [23] . Employing the GIM mechanism results 
in a general expression for the loop integrals which can be written as 

A = V cs V: s [F(x s ) - F(x d )} + V cb V: b [F{x b ) - F{x d )\ , (14) 

with Xi = mf/Mw and F(xd) usually being neglected (except in the 27 case). The s- and 
6-quark contributions are roughly equal as the larger CKM factors compensate for the 
small strange quark mass. The values of the resulting inclusive short distance branching 
fractions, before QCD corrections are applied, are shown in Table 7, along with the 
current experimental bounds.i'0 The corresponding exclusive rates are typically an 
order of magnitude less than the inclusive case. We note that the transition D° — > £ + £~ , 
is helicity suppressed and hence has the smallest branching fraction. The range given for 
this branching fraction, (1 — 20) x 10~ 19 , indicates the effect of varying the parameters 
in the ranges f D = 0.15 - 0.25 GeV and m s = 0.15 - 0.40 GeV. 
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Experimental Limit 

< 1.1 x KT 5 

< 1.3 x 1(T 4 

< 1.0 x 1IT 4 



Decay Mode 



B 



S.D. 



B 



L.D. 



D° 

D° 
D° 



e + e~ 



;i - 20) x 10 



o 



-19 



< 3 x 10 







-15 



D 



o 



77 



IP' 16 
1.4 x 10" 17 



< 3 x 10- 9 



D - 
D° - 
D° - 
D+ 
D+ 



-»• P_ + 7 
-> K*+-f 



< 1.4 x 10~ 4 

< 2.0 x KT 4 



4 x 10~ y 



fewxlO- 10 



< 2 x 10" 5 

< io- 4 

< 2 x 10~ 4 
3 x IO" 7 

1.6 x 10~ 4 



D —> X u + £+£- 

D° -> K°ee/nn 
D° -> p°ee/w 
D + — > 7r + ee//x/x 
L>+ -> K+ee/nn 
D + — > p + nn 



< 1.7 x KT 4 

< 17.0/2.5 x 10~ 4 

< 2.4/4.5 x 10~ 4 

< 250/4.6 x 10~ 5 

< 480/8.5 x IO" 5 

< 5.8 x 10~ 4 



< 2 x 10 

< io- 8 

< 10" 15 



-15 



7T°Z/Z/ 



Z/Z/ 



D + — > 7T + Z/Z/ 



2.0 x 10~ 15 
4.9 x 10~ 16 

4.5 x 10~ 15 
3.9 x 10~ 16 



< 6 x 10~ 16 

< IO" 12 

< 8 x 10~ 16 

< 10" 14 



Table 7. Standard Model predictions for the branching fractions due to short and long distance 
contributions for various rare D meson decays. Also shown are the current experimental limits. 
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The calculation of the long distance branching fractions are plagued with the usual 
hadronic uncertainties and the estimates listed in the table convey an upper limit on the 
size of these effects rather than an actual value. These estimates have been computed 
by considering various intermediate particle states (e.g., n, K, K, r], r]', ttit, or KK) and 
inserting the known rates for the decay of the intermediate particles into the final state 
of interest. In all cases we see that the long distance contributions overwhelm those 
from SM short distance physics 

The radiative decays, D — > X + 7, merit further discussion. One of the goals of a 
high luminosity B physics program is to extract the ratio of CKM elements! V^I/lVy 
from a measurement of B(B — ► p/cu + r y)/B(B — ► K*j). CLEO has placedcJ a bound 
on this ratio of branching fractions of < 0.34, while observing the decay B — > K*^ 
with a branching fraction of (4.5 ± 1.5 ±0.9) x 10~ 5 . This yields a very loose constraint 
on the above ratio of CKM elements. This method of determining the ratio of CKM 
elements depends critically on the assumption that these exclusive decay modes are 
dominated by short distance penguin transitions. If this assumption is false, and the 
long distance contributions to these decays were found to be large, this technique would 
be invalidated. Unfortunately, the theoretical uncertainties associated with computing 
the long distance contributions are sizeable.0 Where separation of the two types of 
contributions is somewhat difficult in the B sector, radiative charm decays provide an 
excellent testing ground. In this case it should be possible to separate out the inclusive 
penguin transitions c — > uj, and determine the rate of the long distance reactions which 
are expected to dominate. For example, the penguin transitions do not contribute to 
D° — ► K* 0/ ~f and it would be a direct measurement of the non-perturbative effects. 
Before QCD corrections are applied, the short distance inclusive rate is very small, 
B(c — » uj) = 1.4 x 10 -17 ; however, the QCD corrections greatly enhance this rate.El 
These corrections are calculated via an operator product expansion, where the effective 
Hamiltonian is evolved at leading logarithmic order from the electroweak scale to down 
to the charm quark scale by the Renormalization Group Equations. This procedure 
mirrors that used for b — > sj, and results in B(c — > uj) = (1.1 — 2.3) x 10~ 5 , where the 
lower(upper) value corresponds to taking the scale of the decay to be 2m c (m c ). (We 
note that these radiative branching fractions have been scaled to semileptonic charm 
decay in order to reduce the CKM and m c uncertainties.) In this case, the rate is 
given almost entirely by operator mixing! The penguin contributions to the exclusive 
channels would then be typically of order 10~ 6 , which is still smaller than the long 
distance estimates in Table 7. Observation of 

7^ tan C (15) 



B(D° -> K*°~f) 



at the few percent level would be a test of the perturbative QCD corrections, or it 
could be a signal for new physics. We note that the predicted values of the branching 
fractions for radiative charm decays are well within reach of the rcF. 

Lepton flavor violating decays, e.g., D° — > / u ± e =F and D — > X + fi ± e T , are strictly 
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forbidden in the SM with massless neutrinos. In a model with massive non-degenerate 
neutrinos and non- vanishing neutrino mixings, such as in four generation models, D° — ► 
/ u ± e T would be mediated by box diagrams with the massive neutrinos being exchanged 
internally. LEP data restricts^ heavy neutrino mixing with e and \i to be l^vef^J 2 < 
7x 10 -6 for a neutrino with mass > 45 GeV. Consistency with this bound constrains 
the branching fraction to be B(D° — > /i e T ) < 6 x 10~ 22 . This same results also 
holds for a heavy singlet neutrino which is not accompanied by a charged lepton. The 
observation of this decay at a larger rate than the above bound would be a clear signal 
for the existence of a different class of models with new physics. 

Examining Table 7, we see that there is a large window of opportunity to discover 
the existence of new physics in rare charm decays. Although the SM short distance con- 
tributions are completely dominated by the long distance effects, there are some modes 
where the size of the two contributions are not that far apart. The observation of any 
of these modes at a larger rate than what is predicted from long distance interactions 
would provide a clear signal for new physics. To demonstrate the magnitude of en- 
hancements that are possible in theories beyond the SM, we consider two examples (i) 
a heavy Q = —1/3 quark contributing to D — > X + £ + £~ and (ii) leptoquark exchange 
mediating D° — > // e^. In the first heavy Q = —1/3 quark may be present, 

e.g., as an iso-doublet fourth generation b' quark, or as a singlet quark in Eq grand 
unified theories. The current boundli on the mass of such an object is my > 85 GeV, 
assuming that it decays via charged current interactions. The heavy quark will then 
participate inside the penguin and box diagrams which mediate the decay and result 
in the branching fractions presented in Fig. 4. The branching fractions are displayed 
as functions of the overall CKM mixing factor for several values of the heavy quark 
mass, and we see that a sizeable enhancement is possible for large values of the mixing. 
A naive estimate in the four generation SM yields! the restrictions \V(&\ < 0.571 and 
\V u b'\ < 0.078. In the second example of new physics contributions we consider lepto- 
quark bosons. Leptoquarks are color triplet particles which couple to a lepton-quark 
pair and are naturally present in many theories beyond the SM which relate leptons 
and quarks at a more fundamental level. We parameterize their a priori couplings as 
Af- „/47r = Fi q a. Leptoquarks can mediate D° — > /i ± e T by tree-level exchange, however 
their contributions are suppressed by angular momentum conservation. From the limit 
B(D° -> /i ± e T ) < lO" 4 , Davidson et alB find 

\J F eu F^ c < 4 x 10" 3 ^ 

where ttllq represents the leptoquark mass. 

Detector studies for the observation of loop-induced decays with leptonic final 
states £ + £~ and X + £ + £~ were performed in Ref . [30] , where it was estimated that the 
level of background was < 1 — 10 events and that a rcF could search for these decays 
with branching fractions down to (3 — 20) x 10~ 8 . This is an improvement of 3 — 5 
orders of magnitude over present limits and puts a rcF within striking range for some 



100 GeV 



(16) 
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of these modes. The radiative D decays will be able to be studied in detail at a rcF. 

6. D° - D° Mixing 

Currently, the best limitsi on D° — D° mixing are from fixed target experiments, 
with xd = Amu/r < 0.083 (where Am^ = irt2 — mi is the mass difference), yielding 
Aitid < 1.3 x 10~ 13 GeV. The bound on the ratio of wrong-sign to right-sign final states 
is r D = T(D° -> £-X)/T(D° -> £+X) < 3.7 x 10" 3 , where 



2 



2 ,AFv 21 

(17) 



in the limit Am D /Y, AT/T < 1. 

The short distance SM contributions to Amp proceed through a box diagram 
with internal d, s, 6-quarks. In this case the external momentum, which is of order m c , 
is communicated to the light quarks in the loop and can not be neglected. The effective 
Hamiltonian is 



n 



Ac=2 
eff 



GfOI 



\V cs V:f {no - mlF 2 O r ) + \V cb V: b \ 2 (l b 3 - mlHO) 



where the Ij represent integra ls§ that are functions of m\ l/Myy and m 2 /m 2 c) and O = 
[«7 M (1 — 7s)c] 2 is the usual mixing operator while O' = [u(l + 7s)c] 2 arises in the 
case of non-vanishing external momentum. The numerical value of the short distance 
contribution is Amp ~ 5 x 10 -18 GeV (taking = 200 MeV). The long distance 
contributions have been computed via two different techniques: (i) the intermediate 
particle dispersive approach (using current data on the intermediate states) yieldsi'S 
Aitid ~ 10~ 4 r ~ 10~ 16 GeV, and (ii) heavy quark effective theory which results!! in 
Am D ~ (1 - 2) x 10~ 5 r ~ 10 _17 GeV. Clearly, the SM predictions lie far below the 
present experimental sensitivity! 

One reason the SM expectations for D° — D° mixing are so small is that there 
are no heavy particles participating in the box diagram to enhance the rate. Hence 
the first extension to the SM that we consider is the addition^ of a heavy Q = —1/3 
quark. We can now neglect the external momentum and Amp is given by the usual 
expression, il 

G 2 F M^m D f2 
6tt 2 

The value of Amp is displayed in this model in Fig. 5a as a function of the overall 
CKM mixing factor for various values of the heavy quark mass. We see that Amp 
approaches the experimental bound for large values of the mixing factor. 

Another simple extension of the SM is to enlarge the Higgs sector by an additional 
doublet as discussed in the above leptonic decay section. First, we examine two-Higgs- 
doublet models which avoid tree-level FCNC by introducing a global symmetry. The 
expression for Amp in these models can be found in Ref. [35]. From the Lagrangian in 



Am, = " F :X fBBn\V^,\ 2 F(m 2 /M^) . (19) 
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Eq. (9) it is clear that Model I will only modify the SM result for Amp, for very small 
values of tan f3, and this region is already excludediE from existing data on b — > S7 
and B® — B d mixing. However, enhancements can occur in Model II for large values of 
tan /?, as demonstrated in Fig. 5b. 

Next we consider the case of extended Higgs sectors without natural flavor con- 
servation. In these models the above requirement of a global symmetry which restricts 
each fermion type to receive mass from only one doublet is replaced^ by approximate 
flavor symmetries which act on the fermion sector. The Yukawa couplings can then 
possess a structure which reflects the observed fermion mass and mixing hierarchy. 
This allows the low-energy FCNC limits to be evaded as the flavor changing couplings 
to the light fermions are small. We employ the Cheng-Sher ansatzpil where the flavor 
changing couplings of the neutral Higgs are Xh°fifj ~ (v^Gf) 1 / 2 '\/WjrfTjA. i j, with the 
rriiCj) being the relevant fermion masses and A^ representing a combination of mixing 
angles. h° can now contribute to Amu through tree-level exchange as well as medi- 
ating D° - D° mixing by h° and t-quark virtual exchange in a box diagram. These 
latter contributions only compete with those from the tree-level process for large val- 
ues of Ajj. In Fig. 6a we show the constraints placed on the parameters of this model 
from the present experimental bound on Atud for both the tree- level and box diagram 
contributions. 

The last contribution to D° — D° mixing that we will discuss here is that of 
scalar leptoquark bosons. They participate in Amp via virtual exchange inside a box 
diagram,c3 together with a charged lepton or neutrino. Assuming that there is no 
leptoquark-GIM mechanism, and taking both exchanged leptons to be the same type, 
we obtain the restriction 

F lc F iu 1967r 2 Am D 
m LQ {A-Ka} D ym D 

where F^ q is defined in the previous section. The resulting bounds in the leptoquark 
coupling-mass plane are presented in Fig. 6b. 

Signatures for D° — D° mixing include like-sign dileptons from D°D° — > i^l^X 
or dual hadronic decays such as D°D° — > K ± n T K ± 7!' T . The hadronic signal can be 
cleanly separated from double Cabbibo suppressed decays at a rcF since quantum 
statistics yield different correlations. This is due to the fact that the double Cabbibo 
suppressed modes can only contribute to the same final states when the D° mesons are 
in a relative S-wave. The interference with the doubly suppressed modes can also be 
use d§ to separate mixing which originates in the mass difference Am d from that which 
arises in the decay Ar^. An extensive study of the double Cabbibo decays will result 
as a product of the search for D° — D° mixing. A thorough Monte Carlo simulation of 
D° — D° mixing at a rcF detector has been performed in Ref . [38] , where the available 
sensitivity for mixing was found to be at the 10~ 5 level. 
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7. CP Violation 



CP violation in the Q = 2/3 quark sector is complimentary to that of the K and 
B systems, but has yet to be explored. The rcF factory could thus provide an important 
first opportunity to explore CP violation in this sector. In the SM, the CKM phase 
is responsible for generating CP violation, and the resulting rates are small. However, 
new sources of CP violating phases could greatly enhance the rates thus rendering CP 
violation in the charm system a sensitive probe for physics beyond the SM. 

CP violation requires the interference of at least two amplitudes with non- vanishing 
relative phases. This can occur indirectly via D° — D° mixing, or directly via asym- 
metries induced in the decay amplitude, or kinematically in final state distributions. 
The first case corresponds to the interference of a D° decaying to a final state / at 
time t, with a D°, which mixes into a D° and then decays to / at time t. This process 
is theoretically clean as the hadronic uncertainties cancel in the asymmetry. However, 
since Amu is extremely small in the SM the induced CP violation is negligible. If new 
physics were to enhance D° — D° mixing, as seen to occur in the previous section for 
some models, then this mechanism could yield sizeable CP violating effects. 

Direct CP violation in charm meson decays yields slightly more promising results. 
Here, only decay amplitudes with two separate weak phases and two different strong 
phases will contribute. This can be easily seen as follows. Let us assume that the decay 
amplitude to final state / has the form 

A f = A ie iSl + A 2 e i52 , (21) 

with Ax,2 being the two amplitudes after the strong phases 8\ t2 have been factored 
out. For the CP conjugate amplitude, the weak phases are conjugated, but the strong 
phases are not. The CP asymmetry is then given by 



\A f \ 


2 _ 


\*f\ 


2 




21hi{A* 1 A 2 )sm(8 1 -8 2 ) 


\A f \ 


2 + \ 




2 


Ail 


2 + \ 


A 2 \ 


2 + 2-Re{A* 1 A 2 )cos{5 l -5 2 ) 



(22) 



which clearly vanishes if A\ )2 contain the same weak phase and if 8\ = 8 2 . Before 
estimating the typical size of this asymmetry in the SM, we first note that in contrast 
to B decays, the branching fractions for the relevant modes, i.e., 7r + 7r~, K + K~, etc., 
are rather sizeable in the charm system, and for once, the large effects of final state 
interactions are welcomed! Of course, this makes numerical predictions difficult due to 
the hadronic uncertainties associated with final state interaction phases. All tree-level 
interactions in D decays contain the same weak phase, hence they must interfere with 
penguin mediated processes in order to obtain the requisite two independent weak 
phases. Since Cabibbo favored modes are not mediated via penguin type interactions, 
only Cabibbo suppressed decays exhibit direct CP violation. Although the relative size 
of the tree and penguin contributions has not yet been determined, an upper bound 
for the size of the CP asymmetry is estimated!!®!! to be acp < 10~ 3 . However, there 
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could be a large enhancement from the strong interactions (which could be provided 
by, e.g., nearby resonances) and typical results in this case are acp ~ few x 1CT 3 . 

It is possible to obtain two separate weak phases in tree-level amplitudes in D s 
decays, i.e., via interference between the spectator and annihilation diagrams. If the 
annihilation processes are not suppressed relative to the spectator case, then CP asym- 
metries of order 10~ 3 are feasible. 

Kinematic CP violation signals could occur, for example, in the decays D — > VV , 
which are described by more than one amplitude. Here, it is possible to construct CP- 
odd triple product correlations between the two polarizations and one of the momenta, 
e.g., (k • ei x 62). Since final state interactions can induce a non-zero value of the triple 
product correlation and hence mask CP violation, one must evaluate 

Njk ■ ex x 6 2 > 0) - Njk ■ ex x 6 2 < 0) 

7V / = ^ > K 26 ) 

IV total 

as well as the corresponding quantity for Nf, which vanish if CP is conserved. If Nf 
and Nj are determined to be non-zero then there is a clear signal of CP violation. 
Similar' correlations have also been diseussedB for semileptonic decays. 

An interesting example of the potential size of CP violating effects from new 
physics is that of left-right symmetric models.@@ In this case reasonably large values 
for CP asymmetries can be obtained for the Cabbibo allowed decay modes. This occurs 
due to the existence of an additional amplitude from the Wr exchange, which carries 
a different weak phase from that of the Wl mediated decay. The estimated values 
of the CP asymmetries in these models is of order 0.01. Compared to the vanishing 
asymmetry obtained in the SM, this would provide a sizeable and clear signature for 
new physics. 

The experimental feasibility of discovering CP violation at a rcF has been exam- 
ined in detail as in the case of D° — D° mixing, where Monte Carlo results indicate 
that the backgrounds to the two-body hadronic decays are manageable. These authors 
find that a rcF should be able to probe direct CP violating asymmetries in D decay 
at the level of 10~ 3 , which just touches the range predicted within the SM. 

8. Summary 

In summary we see that the charm physics program at a Tau-Charm factory is 
diverse and robust. We have seen that (i) it is possible to greatly improve the precision 
over present data samples on the charm system. Most decay modes as well as the 
pseudoscalar decay constants and the CKM elements governing charm-quark decay 
can be measured at a level of order 1%. This would greatly enhance our understanding 
of the underlying QCD interactions which govern these processes. In some decay modes 
and other interactions the rcF could provide the first opportunity for observation, (ii) 
Loop-induced reactions could be probed at a substantially increased level of sensitivity 
allowing for a wide window of opportunity to discover new physics, (iii) CP violation 
could finally be explored in the up-quark sector, where the capabilities of a rcF in 
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one year of running at design luminosity just touches the range of SM predictions. 
Several years of running (or increased luminosity) could study CP violation at the SM 
level. This wide physics potential and the high precision which can be achieved would 
place measurements in the charm sector at an equal level with those in the down-quark 
systems. 

The rcF has the potential to address fundamental questions of the SM in a 
complimentary fashion to the on-going extensive studies of the down-quark sector. 
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Fig. 1: Distribution^ of difference between missing energy and momentum for Cabibbo 
allowed and suppressed semileptonic decays. The shaded regions correspond to the 
background levels assuming (a-b) lead-proportional tube calorimetry and (c-d) Cesium 
iodide crystal calorimetry. 

Fig. 2: Missing mass distributions!!! for (a) D + — > /U + z/, (b) Df — > (c) D + — > t + u 
with r — > fiuu, and (d) Df — > t + v with r — > ez/z/. The shaded regions correspond to 
the background levels. 

Fig. 3: Reconstruction! of the decay D s — > ?77r + 7r + 7r _ with 77 — > 7r + 7r _ 7r and 7r° — > 77. 

Fig. 4: Branching fraction for D — > X + in the four generation SM as a function of 
the CKM mixing factor, with the solid, dashed, dotted, dash-dotted curve correspond- 
ing to m v = 100,200,300,400 GeV, respectively. 

Fig. 5: Amn in (a) the four generation SM with the same labeling as in Fig. 4, (b) 
in two-Higgs-doublet model II as a function of tan f3 with, from top to bottom, the solid, 
dashed, dotted, dash-dotted, solid curve representing m H ± = 50, 100, 250, 500, 1000 GeV. 
The solid horizontal line corresponds to the present experimental limit. 

Fig. 6: (a) Constraints in the mass-mixing factor plane from Am^ in the flavor changing 
Higgs model described in the text due to the tree-level process (solid curve) and the box 
diagram (dashed), (b) Constraints in the leptoquark coupling- mass plane from Am^. 
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